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After Sakharov's seminal paper on the conditions for obtaining a matter asymmetric 
universe several authors have attempted to realise these conditions in the context of different 
models[|l|. One approach is to create a lepton asymmetric universe by the asymmetric decay 
of heavy scalars or Majorana neutrinos and to convert the lepton asymmetry into a baryon 
asymmetry by B+L violating sphaleron processes. In the original papers on leptogenesis CP 
violation was realised by the interference between the tree level diagram and the one loop 
vertex correction[§. Later it was realised that there was an additional contribution from the 
one loop self-energy correction to the external heavy scalar or neutrino0, |]. 

The self-energy contribution has been calculated differently by different authors. Refs. |@, 
^ obtain the asymmetry by considering the interference of the tree level diagram for decay 
and the diagram containing the one loop self-energy correction. However this approach 
breaks down in the degenerate heavy particle mass limit. In refs.[|], ^ ||, ||, [10| the authors 
obtain the one loop effective Hamiltonian, or more precisely, the one loop effective mass 
matrix for the heavy species. Because particles in the self-energy loop go on shell the effective 
mass matrix is non- Hermit ian. This, coupled with CP violation, results in different evolution 
of the scalars and their antiparticles, or of the right-handed and left-handed neutrinos, and 
ultimately gives rise to an asymmetry when the heavy particles decay. In refs.[^ |T0| the 
authors find an enhancement in the degenerate scalar or neutrino mass limit |Tl|] . The authors 
of refs. 0, 0, H, 1^ and |I0| differ in their use of the effective Hamiltonian. In this paper we 



compare and comment on both approaches and on the approach in refs. p^, |T3| where the 
authors use resummed propagators to evaluate cross sections and decay rates. We shall 
discuss below models involving decays of heavy Majorana neutrinos. However our comments 
also apply to models with heavy scalars. In refs.p, |10| the authors consider models with 
charged scalar fields; to facilitate discussion we have modified their definitions to suit the 
Majorana fields Na- 

In section 1 we discuss some subtleties associated with a non-Hermitian Hamiltonian, 
particularly with respect to the normalisation of eigenstates of the Hamiltonian. In section 2 



we compare the approaches of refs.0 and ||T^ both of whom use the effective Hamiltonian. 
In section 3 we discuss some field theoretic concerns regarding the definition of the decay 
amplitude for eigenkets of a non-Hermitian Hamiltonian and the existence of a Fock space 
for mixed fields. In section 4 we compare the matrices used to diagonalise the effective 
Hamiltonian or projections of the resummed propagator in refs.[^ ^ We also comment 
on the approach of ref.|T^. We re-evaluate the asymmetry in the small mass difference limit. 
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In section 5 we summarize our results. 



1 Diagonalisation of the effective mass matrix 



Let us start with the lagrangian 



int 



J2M,[{NmYNm + NmiNR.,y] 

i 



Ri) 



(1^ 



as in ref. 0]. A^f^ = {NmY and A^'^ = CN . 11 are hght leptons and is a scalar field. For 
concreteness the /£ can be the left-handed lepton doublet of the Standard Model, can be 
the Standard Model Higgs doublet. We supress the SU(2) indices carried by the lepton and 
the Higgs fields which are contracted by the eab tensor, a and i are generation indices and 
we assume one heavy Majorana field per generation of light leptons. For simplicity we shall 
work with two generations. The lagrangian may also be rewritten in terms of rji = Nm + N^j^. 
In what follows we shall replace N^j^ by Nn. 

Since the right handed neutrinos Nm decay to leptons while the left-handed neutrinos Nn 
decay to antileptons, the lepton number violation is in the mass term for the heavy neutrinos. 
Note that the last two lines in the lagrangian above are equivalent as X1X2 = XiXi- O^ig 
loop effects give rise to mixing between heavy neutrino flavours as shown in flg. 1. The one 
loop effective mass term due to the diagrams in flg. 1 is given by 
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Ignoring the masses of the light leptons and of the scalar fleld, the expressions for ifj^ and 
ifjj are given in ref.|0] as 



a a 



{9, 



dis _n'^^^} 
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(3) 
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Figure 1: One loop contributions to the mass matrix 



(9^ - (5) 



and 



where = ^1^. 

The effective mass matrix Mg// = Mq + if in eq. (j^) can be written as M — |r where the 
dispersive part of H (i.e. terms proportional to is included in M and the absorptive part 
(proportional to (yf^^j) is included in V. The absorptive part reflects the fact that the particles 
in the loop have gone on shell. Both M and F are Hermitian but not M^jf. This is analogous 
to the formulation used to obtain an effective mass matrix for neutral kaons|n]. The one 



loop corrections in M can be absorbed away by renormalisation. (The off-diagonal elements 
of M are renormalised by non-diagonal elements of the wave function counterter ms ||1 5|] . ) 
However, there are no absorptive counterterms ||1 5|| , nor are they needed for the absorptive 
parts are finite. It will be the mixing via the absorptive terms that will be important for 
creating the lepton asymmetry. 

Since the effective mass matrix is non-Hermitian it can not be diagonalised by a unitary 
matrix, i.e. VM^ffV^^ = M^, where is diagonal, but ^ V'^ . The columns of 
V^^ are the eigenkets of Me// and the rows of V are the suitably normalised eigenbras of 
Meff, i.e., xMeff = Xx. (The subtleties of normalisation are discussed below.) If Me// were 
Hermitian the eigenbras of Me// would be the Hermitian conjugates of the eigenkets and 
would equal V"^ . This is not true in our case, though the bra-eigenvalues and ket-eigenvalues 
are the same. 

Therefore we rewrite the mass term of eq. (0) as 

NMeffN = i^Mo^jj = i>M'i} (6) 

where if) = VN and = NV~'^, and only M^ is diagonal above. Besides the choice 
mentioned above one can diagonalise Me// by a bi-unitary transformation, i.e, U^Mf,ffV = 
Md, where U and V are unitary. However in this case the kinetic energy term will not be 
diagonal in the new (■?/', i') basis. 

If one can expand the ipjip and N fields in terms of annihilation and creation operators 
(later we comment on this assumption), then {ipc) = \Na)V^^, where {ipc) represents a particle 
created by the action of the creation operator in ipcix) = Na{x)V^^ on the vacuum and \Na) 
represents a particle created by the creation operator in Na{x) acting on the vacuum. Note 
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that while Ea \Na){Na\ = 1, Ea \A){^a\ ^ 1 wheie (^,1 = l^aV []• Instead Ea \i'a){^a\ = 1 
and {ipali^a) = 1 where {ipa\ represents a particle created by the action of the destruction 
operator in ip{x) acting to the left on the vacuum^. 



2 Obtaining the lepton asymmetry 



Refs. 1^, 1^, 1^, 3 and [|T^ differ in their choice of the initial state relevant for the calculation of 
the lepton asymmetry. Let e be a measure of CP violation due to the self-energy correction. 
In refs.|P, 0, ||, || the authors define 



(7) 



while the definition of the asymmetry parameter in ref . |TD 



IS 



r 



r 



+ r 



(8) 



where, as stated earlier, the definitions of the asymmetry parameter in refs. fg, |I0[ have been 
modified for a scenario involving Majorana neutrinos. Below we compare the approaches of 
refs.0 and 0. 

Before we proceed let us reiterate here the mechanism of the creation of the lepton asym- 
metry. Starting with two relativistic neutrino species the second (heavier) species becomes 
non-relativistic at tnr2- After t„r2 the physical states are the eigenstates of an effective 
Hamiltonian for one relativistic neutrino species with a Lorentz suppressed mass and one 
non-relativistic neutrino species. They can be represented by A'^i and ip2 fields, where ip2 is a 
linear combination of the Na (but is not tp2)- We assume the heavier decay while the Ni 
neutrinos are in equilibrium. At time teq the lighter neutrinos A^^i go out of equilibrium and 
later at tnr they become non-relativistic, i.e., we assume K = Mi /Teq <^ 1, where Teq is the 
temperature at teq- Thus at tnr one will have a universe with only A^^i neutrinos with equal 
numbers of right-handed and left-handed neutrinos, Nm and A''^;^. After tnr the effective 



^In ref. ^ the authors discuss the fact that because of the non-Hermitian nature of the effective mass 
matrix -0 ^ i/j'^. They state that ip still represents a Majorana field and hence we sum only over particle 
states. In models involving charged scalars one certainly needs to consider both particle and antiparticle 
states. In that case the sum over all states will include particle and antiparticle states. 

^ The normalisation of ip states has also been discussed elsewhere, as for example, in ref. ||l6| and references 
thereof. 
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Hamiltonian is as in eq. and the ipa are the eigenstates of the effective Hamiltonian. 
The state of the universe with equal numbers of Nm and Nil can now be described as a 
hnear combination of the ipa- More precisely, a state \Nib) (or \Nil)) at tnr can be rewritten 
in terms of \ipa)- As in the standard neutrino flavour oscillation scenario, the right- and 
left-handed content of an originally |A^i_r) or IA'^il) state will now start to oscillate. One 
can see from eq. (|) that the amplitude for \Nir) |A^2l) is proportional to H12 while the 
amplitude for \Nil) — * |A^2_r) is proportional to Hi2. (The amplitudes for \Niji^l) iNn^ji) 
are equal and proportional to Mi + Hu.) Because of CP violation H12 7^ H12. Therefore, 
though one starts with an equal number of Niji and Nn, more Nn may transform into 
N2R than Nm into N2L and one ends up with an universe with an unequal content of Nm 
plus A^2_R and A^i^ plus A^2L- When the neutrinos subsequently decay one obtains a lepton 
number asymmetry since right-handed neutrinos decay to leptons and left-handed neutrinos 
decay to anti- leptons. (We ignore here asymmetry due to the vertex correction.) 

The formalism of ref.|jlO| starts with a universe with equal numbers of Nifi and Nn at 
tnr, or, equivalently, with states |A^ir) and \Nii,) with equal probabihties of existing at tnr- 
The asymmetry obtained reflects the oscillations in the right-left content that takes place 
prior to decay. On the other hand, the final asymmetry is obtained in ref.|^ by considering 
only {ipi). One may argue that typically one does consider the asymmetry generated by only 
the lightest neutrino. But what that entails in this scenario is that the heavier ip2 neutrino 
decays while the interactions of the lighter A^^i are still in equilibrium, leading to the erasure 
of any asymmetry generated in the decay of tp2. Subsequently the lightest neutrino A^^i goes 
out of equilibrium, becomes non-relativistic and is re-expressed as a linear combination of 
both \ipa) states. Thus one has to consider the lightest Na species and not the lightest ipa 
species as the initial state at tnr- Nevertheless, since {ipi) is an asymmetric mixture of \NaR) 
and {Nai) its decay gives a lepton asymmetry. 



Ref. also argues against starting with the ipa states as these are "propagation" states 
which the authors claim are not physical. However our observations are based on a cos- 
mological argument rather than a field theoretic argument. (We shall also state the field 



theoretic argument later.) The importance of the initial state is also discussed in ref. [|T0 



Our arguments above do not apply to scenarios where K ^ 1, i.e., when the neutrinos 
are non-relativistic when they go out of equilibrium. In such cases, from a cosmological 
standpoint, e^^ will be the relevant asymmetry parameter. 

In the limit of large mass splittings, i.e., \Mi — M2I ^ the asymmetry obtained 



7 



in ref.[0 is similar to that obtained in ref.|0]. We find that the approximations needed to 
obtain this result, i.e., ignoring terms proportional to W1W2, W2W1 and 11^2^ in eq. (17) 
of ref. and assuming a small mixing angle is equivalent to ignoring the asymmetry from 
the interference between ipi and ip2 and from 1IJ2 decay and to setting A^^i ^ ipi- However in 
the small mass splitting case the asymmetry obtained in refs.[^ |T^ are not similar, as we 
show later. 

Our arguments imply that one has a mixed state (as opposed to a pure state) of \Nifj) and 
\Nil) with equal probabilities at teq. Thereafter one evolves the \Nir^l) states independently. 
Since the mass term is not significant when A^^i is relativistic there will be no significant mixing 
till tnr- After tnr both INiji^l) can be rewritten as pure states of and \ip2) and hence one 
gets interference between {ipi) and \ip2) states which is necessary for left-right oscillations. 
However if one argues that one instead has a mixed state of lipa) at tnr then one would 
evaluate e^^ as in ref.0. In that case the initial number densities of the ipa particles that 
enter into the solution of the Botzmann equation would be weighted by the probabilities of 
obtaining the ipa at tnr, given by KipalNi)]"^. One would then include the contributions of 
both ipa to the total asymmetry since the asymmetry generated by the decay of ip2 will not 
be erased as ipi interactions are out of equilibrium at tnr- 



3 Some field theoretic concerns 



We now raise more fundamental issues with both the approaches using the effective Hamil- 
tonian. The asymmetry parameter is defined in refs. |^, |^, |^, ^ using decay rates for rotated 
fields ipa with complex masses. It is shown in ref. (in the context of kaon physics) that 
the lack of orthonormality of the \tpa) states raises ambiguities in the normalisation of states 
making it difficult to interpret {flip a) as the probabihty amplitude for ipa f- If one nor- 
malises l-ipa) as |V'a)/'\/ (V^alV^a) then the state is scale dependent. This is because \ijja) — > c|'?/'a) 
implies {ipa\ —>■ ("^Aalc"^ since {ip) = \N)V~^ and = V{N\. On the other hand, if one 
normalises the states as \i)a) / \J {i'a\i'a) then E/ K/|V^a)P/[(^a|^a)](/|/) 7^ 1- This can 
be better understood by considering the quantum mechanical example of kaon mixing with 
iV'a) corresponding to \Ki) (or \Ks)) and (/| corresponding to (-R'l.s'I and {Kl^s\- A similar 
argument can be extended to field theoretic scenarios. If ipa are intermediate particles, i.e., 
if one includes a production amplitude for the ipa-, the above concerns do not apply If 
one does include a production amplitude for the ipa in e^^ it will cancel out in the numerator 
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and denominator leaving behind the asymmetry parameter as defined in refs.[^ ^ ^ ^. 



In ref. |10] the asymmetry parameter is defined in terms of decay rates of the mixed Na 
fields. It has been pointed out in ref. that in general there does not exist a Fock space for 
mixed fields. If one assumes that the annihilation operators of the ipa fields and the creation 
operators of the ipa fields obey canonical anticommutation relations then the annihilation 
and creation operators of the Na fields satisfy, for example, |jTB| 

E'^ - h\k\ 
'2EL 



{A{k, h, t), A]{k', h', t)} = S'{k - k')5w E Vi-a'Va, \J ' (9) 



where h, h' = ±1 represent helicity, E'^ = (|kp + — imaVa)^ , the sum is over the mass 
eigenstates, and A^j include spinor elements as in ref.|T8| and we have used = V'^ip 
and N = ipV. These anticommutation relations are non-diagonal in the i, j indices and are 
hence not canonical. Therefore if iV = ipV, \N) = \^p)V is not a well defined statement in 
general. Furthermore, the two cases in which one can define an approximate Fock space for 
the N fields, namely, the extreme relativistic limit and almost degenerate masses with real 
V [|T^ do not apply to our problem. 

In the limit that all the neutrino species are non-relativistic, as is the case for t > tnr, 
the above commutation relations are diagonal in the i,j indices. However, as pointed out in 



ref. [|18| , the annihilation and creation operators of the N fields are no longer time- independent 
in this case. But if one starts with an Na state as the initial state at tnr and evolves it as a 
linear combination of ipa states then this may not be a problem since one requires that the 



relation |A^) = \^p)V hold at only one instant, namely, at t^r- Note that in ref.[jTO[ the states 
and |A^) are in co-ordinate space. The creation and annihilation operators for the mixed 
fields in co-ordinate space, namely, N{x) and N{x), do satisfy diagonal anticommutation 
relations. 



4 Final comments 

It was pointed out by Veltman that it is inappropriate to have unstable particles as asymp- 
totic in or out statesj^. To avoid this problem different effective approaches have been 



adopted to obtain the one loop decay amplitude. In refs. |^, |, ^, |T^ the authors use an 
one loop effective Hamiltonian. In ref. the decay amplitude is obtained by splitting a 
two particle scattering amplitude with the unstable particle as the intermediate particle. In 
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ref.[|T3| one obtains the decay amplitude from an effective LSZ formulation. We discuss here 
the differences and similarities in the above approaches to obtain the asymmetry.^ 

Broadly speaking, one can further divide the different approaches into those that use 
unrotated fields Na to define the initial state, such as refs.[|T^, [T^ and those that work with 
rotated fields to define the initial state, such as refs.]^, |^ |^, ^ 00- Though refs.[|^, ^, 
deal with rotated fields their approaches are somewhat different. Below we first compare the 
rotating matrices for all these approaches in the small and large mass difference limits. We 



also comment on ref . |I3| . We finally re-evaluate the asymmetry in the small mass difference 
limit. 

In ref. the authors obtain the eigenvalues and eigenvectors of the effective mass matrix 
in eq. (Q). In ref. the authors obtain the eigenvalues and eigenkets of the two off-diagonal 
blocks of Me//. In general, the eigenvalues of the two off-diagonal blocks of the matrix do 
not equal each other or the eigenvalues of the full matrix and the eigenkets of the off-diagonal 
blocks describe particles different from those related to the eigenkets of the full matrix. 

One can show that the eigenvalues of the two off-diagonal blocks are the same as that 
of the larger matrix, and the eigenvectors of the larger matrix are trivial combinations of 
corresponding eigenvectors of the smaller blocks when the two off-diagonal blocks are the 
same. In the limit of a small mass difference |Mi — M2I ^ |-f^i2|, H12 ~ H12 and we find that 
the diagonalising matrices are equivalent and one gets the same asymmetry. For example, 
the matrix relating the right handed field projections Napt, and ipaR is given by 



1 _i I 



(10) 



for both approaches. The above matrix is obtained from V'-^ of Section 1 and we have not 
shown the normalisation. M- = Mi + Ha. However in the limit of a large mass difference 
I Ml + Hii — M2 — H22\ ^ |-f^i2| the relation between NaR and ipaR is given by 



/ -, _M[Hi2+M^Hi2 



(11) 



•^In section 7.3 of ref. [g2j a further approach to obtain the decay rate for unstable particles using the 
optical theorem is given. In ref. 1 23 a modified leptogenesis scenario in which the Majorana neutrinos are 
virtual intermediate particles is considered to get around some of the issues discussed here. 

^It is pointed out in section 5 of ref.|^ that in some cases diagonalising the effective mass matrix to 
obtain the ipa fields may not be possible. 



10 



and 

1 gl2 
Hl2 I 



(12) 



for refs.[0] and |^ respectively. One can see that they are not equivalent. We find that 
working with the 2x2 matrices in this limit gives an asymmetry that is approximately 1/2 of 
the asymmetry obtained using the 4x4 approach. This can be traced to the absence of the 
contribution of the term proportional to M[Hi2. Once again, if hai are real H12 = H12 and 
the diagonalising matrices of the two approaches become identical. 

In ref.Q the authors also consider the eigenvalues of the two blocks of a 4x4 effective 
mass matrix. However our above arguments do not apply to this case as the two blocks are 
along the diagonal and are transposes of each other, necessitating that the eigenvalues of 
the two 2x2 blocks are the eigenvalues of the larger matrix and that the eigenvalues of each 
block are the same, and the eigenvectors of the larger matrix are a trivial combination of 
the eigenvectors of the blocks. 

Ref. 1^ diagonalises the mass matrix and obtains the relationship between \ipa) and \Na) 
states given by {ip) = \N)V~^. On the other hand, ref. |]12| diagonalises projections of the 
resummed propagator, Srji = {NjiNl) and Sll = (^l^r)- Using this to factorise the 
amplitude for l*(j) and l*(f) ip 1(f)* they obtain the amplitudes for l*(j) 

in the large mass difference limitQ. In the limit of a large mass difference, \M[ — M2I ^ \Hi2\-, 
the matrix relating NaR and ipaR fields for ref. |jl2| is given by 



/ 1 _ 2M2^i2 

M'(-M'^ 

2M1H12 1 



(13) 



v 



M(-M^ 



Though the matrix above differs from that for ref. the ultimate asymmetry is the same as 

that of ref.0 because the factor of 2 in the off-diagonal entries above compensates for the 

presence of the terms proportional to H12 in eq. ( pH] ) We have also obtained the matrix 

relating NaL,R and ipaL,R fields in the limit of a small mass difference, \M[ — M2I ^ \H12\1 

by diagonalising the projections of the resummed propagator as in ref.|]T2|]. We obtain the 

same matrix as refs.|]^, ^, as given in eq. (pi]|)f|. 

^In section 9 of ref.|2l| a similar factorisation of two particle scattering amplitudes was done (in a different 
context) to obtain the amplitude for decays of unrotated fields. 

^Note that the lagrangian of ref. as given in eq. (|^) is twice that of ref. . 

^In general Srr and Sll do not have the same poles making it non-trivial to define the particle mass. 
We do not address such field theoretic subtleties. We only mention that in the limit that \M[ — Mjl ^ \Hi2\ 
Srr and Sll do have the same pole. 
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In ref . the author obtains the resummed propagator as in ref . [0 but does not diago- 
nahse it and obtains the amphtudes for A^i — ^ l*(j) by by truncating the Green's function 
in momentum space with (S'n)^"'^ = fi — Mi + Sn — t^i2[i> — M2 + S22]^"'^S2i, where the 
self-energy corrections to the propagator Sj^ = iE^j** are defined in eq. (4.2) of ref. |T^. In 



the LSZ formalism the truncating operator defines the corresponding incoming or outgoing 
particle. This implies that the incoming neutrino has a mass that is the (complex) pole of 
the component Su and not of the full resummed propagator S, i.e., one is including only one 
loop diagonal corrections to define the neutrino mass and excluding off-diagonal corrections. 
Excluding mixing in the incoming state allows one to use the classical Boltzmann equation 



as discussed later. The derivation of the LSZ formula in momentum space used in ref. [|13 1 
involves integrations over time in intermediate steps, whose range of integration must be 
modified from (—00, +00) to (0, +00) if the neutrino fields have a complex mas^. However 
the decay amplitude obtained in ref. can be embedded in an amplitude involving only 
stable particles in asymptotic states (analogous to the approach in section 9 of ref. |^), if 



one equates = — Mi + Su — f^ui^ — M2 + S22) ""^^21] with uiUi/{...) and again ui 
and Ui satisfy only a diagonal Dirac equation]^. ("..." is approximately equal to iMiTi where 
Fi includes width effects close to Ni production.) 

We reiterate here the nature of the incoming particle in various approaches. In ref.|Q 
the incoming particle is an unmixed particle corresponding to the original lagrangian of the 
theory, i.e., Nil^r of eq. (|1|). In ref.[|I^ the incoming particle is a particle corresponding to 
the mixed field of the one loop effective theory, i.e, Ni^^ji of eq. (@). In ref.H, 0, H, H the 
incoming particle is a particle corresponding to the rotated field of the one loop effective 
theory, i.e., ipi of eq. (§) or its equivalent. In ref. the incoming particle corresponds to 
the rotated field obtained by diagonalising projections of the resummed propagator. The 



incoming particle has a one loop effective mass. In ref. |]13| the incoming particle is unrotated 

but has a one loop (complex) effective mass corresponding to the pole of the Su component of 

the resummed propagator. In the limit that mixing is small compared to the mass difference, 

i.e, 1-^12 1 ^ \M[ — M2I the amplitudes for decay in refs.p, [12|, |13[ are similar though the 

spinor u in the amplitude satisfies uu = ^ + m, or uu = ^ + m, with a different m for each 

approach. As implied earlier the amplitude for ref.|^ differs in form from the others but 

^Integrating over a semi-infinite time interval has been taken into account in ref. [lC|] (and, for example, 
in ref.[|^ in a different context). 

®For a field with a complex mass m, uu, and not uu, equals ^ + m. 



12 



refs. 1^, 0, [T^ give approximately the same asymmetry in this hmit[^ 



In the opposite limit that \M[ — M2I ^ |-f^i2| the rotating matrices of refs.J^, 0, 112 



are equivalent and give the same asymmetry. It was first shown in ref. [0 that there is an 
enhancement in the asymmetry when |Mi — M2I < |-f^i2| and \Hii—H22\ <^ |-f^i2 + -f^i2|- We 
have reevaluated the asymmetry in this limit and obtained 



Im{hl^K2hlihp2) 



in V^[2\hs2\^ + Reih;,h^2)] + A{gf4YRe\K^K2)\hsi + hs2^^ 



(14) 



where 77 = (M2 — Mi) /Mi and the sum over a, /3, 7, 5 and e is implied]^. The expression for 
the asymmetry obtained in ref.|T3| in the limit |Mf — M|| <^ Mf is 

_ Im{h*^iha2h}ihf32) tn , . 

where = {Mf - M|)/(MiM2) ^ -2r] and Aij = Ea haih^^j / {16n) . In obtaining the 
asymmetry above, in the invariant amphtude was substituted by Mi. To be consistent with 
the truncation of the Greens's function should be replaced by the pole of Su. It is not trivial 
to solve for the pole. Using a naive substitution of by Mi — Su + Si2[Mi — M2 + S22]~"'^S2i, 
where are evaluated at |^ = Mi and ignoring terms of order h^^ at the amplitude level, 



as in ref. |]r3|, one obtainsQ 

, _ Im{hliK2h}ihp2) rjv(l -g) 

""^^ 87r|/i5i|2 r^(l -a)2 + 4(^22- An + aA22)2' ^ ' 

where a = {A^l'^Mf /[{Mi - M2f + Al^Ml]. In the case of hierarchical Yukawa couplings 
the above result reduces to the asymmetry in eq. (pIS]). In the limit that An ~ A22, or 
l^ail ~ \ha2\1 and 7] ~ ^22, which is comparable to the limit in which the enhancement was 
seen in ref.[0, a ~ 1/2. Then the asymmetry reduces to 

, _ Im{h*^^K2hlihp2) 2rN , . 

»7r|ft5i|^ r^+4/l22 



^"^Whcn the incoming particle has a complex mass there are complications in obtaining the decay rate 
by squaring the decay amplitude because uu =^ ;^ + m. However this does not affect the asymmetry which 
involves a ratio of decay rates. 

^^This result differs from that of ref.0 but we have been informed by Dr. U. Sarkar that there is a 
correction to the calculaton of e^^ in ref. 0| in this limit. 

^^Finite temperature effects, which we have not included, can induce changes in particle masses which can 
be relevant lE5|. 
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One can see that the asymmetry in eq. is different from that in eq. (p!?]). This is not 
surprising as eq. ( p!^ represents the asymmetry in the decay of a. ip particle while eq. ( p!7D 
reflects the asymmetry in the decay of a particle represented by an unrotated field with a 
complex mass. In the limit of large mixing these need not be equivalent. 

We have also obtained the asymmetry in the small mass difference limit for the model 
involving sneutrino decays in ref.|T^. In the limit that |Mf — M|| -C and |r^]^ — rl2\ ^ 
|r^2l (i^ the notation of ref.|jTO|), we find 

^ Im{hl^K2h}ihp2) rjv , . 

'''' T^{\hsi\^ + \hs2?) Tl + lQ[{A22-AuY+Um^A,2V ^ ' 

where Aij and are defined as above. While comparing this result with the asymmetry 
obtained in other works one must keep in mind that eq. (|T^) is the asymmetry obtained 
from sneutrino decays and that the model in ref.[|I^ has twice as many decay channels for 
the sneutrinos as compared to the neutrinos in other works. 

We would also like to add a note of caution when one uses the effective Hamiltonian 
approach and the Boltzmann equation to obtain the lepton asymmetry. The Boltzmann 
equation which is used to calculate the number densities of various species assumes that 
the nature of the species does not change over time. The l.h.s. of the Boltzmann equation 
is the time derivative of the number density of a species and one has to be careful as to 
exactly what one means by nAr(t) or n^{t) at times when A^^ or ipa may not be the physical 
states. This comment applies to scenarios where one starts with states at tnr- which evolve 



subsequently as ip states, as well as to ref . |T9[ where one obtains the asymmetry by evolving 
the %l) number density from when the neutrinos are relativistic. 



It has been argued in ref. [^] that decoherence effects due to interactions with the thermal 
universe justify excluding mixing effects in the incoming neutrino state. Excluding mixing 
in the incoming state allows one to work with a classical (diagonal) Boltzmann equation. 



5 Summary 

In conclusion, we have argued that in leptogenesis scenarios where the heavy neutrinos or 
scalars are relativistic when they go out of equilibrium one should calculate the asymmetry in 
the decay of the lightest unrotated species. We have pointed out certain issues related to the 
definition of the decay amplitude for eigenkets of the non-Hermitian effective Hamiltonian 
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and the existence of a Fock space for mixed fields. As different authors have used various 
approaches to rotate the mixed fields we have compared the rotating matrices of refs. ^ ^, |12[ 
and discussed their similarities and differences. Finally we have re-evaluated the asymmetry 
in the small mass difference limit using different approaches. 
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